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ABSTRACT 

We present a 40GHz (7.5 mm) raster scan image of a 3.6° x 2° region centered on the low redshift 
(z = 0.055) cluster of galaxies Abell 3667. The cluster was observed during the Antarctic winter of 
1999 using the Corona instrument (15.7' FWHM beam) on the Viper Telescope at the South Pole. The 
Corona image of A3667 is one of the first direct (i.e. rather than interferometer) thermal Sunyaev- 
Zel'dovich effect images of a low redshift cluster. The brightness temperature decrement at the X-ray 
centroid (20H2 m 28.9 s , -56°49'51" J2000) was measured to be AT C mb = -154fiK. We have used the 
40GHz map of A3667 in conjunction with a deep ROSAT PSPC (X-ray) image of the cluster, to make 
a measurement of the Hubble Constant. We find Hq = 64^30 km s _1 Mpc^ 1 (68% confidence interval). 
Our Hq calculation assumes that the cluster can be described using an isothermal, tri-axial ellipsoidal, 
/3-model and includes several new analysis techniques including an automated method to remove point 
sources from X-ray images with variable point spread functions, and an efficient method for determining 
the errors in multi-parameter maximum likelihood analyzes. The large errors on the Hq measurement 
are primarily due to the statistical noise in the Corona image. We plan to increase the precision of our 
measurement by including additional clusters in our analysis and by increasing the sensitivity of the 
Viper SZE maps. 



1. INTRODUCTION 

The Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 
1972, SZE hereafter) describes the inverse Compton scat- 
tering of cosmic microwave background (CMB) photons by 
energetic free electrons. For instance, the random thermal 
motions of electrons trapped in the potential wells of clus- 
ters of galaxies result in a frequency dependent change in 
CMB intensity known as the thermal SZE. Similarly, the 
bulk peculiar motion of these electrons results in frequency 
independent change known as the kinetic SZE. Observa- 
tions of the thermal and kinetic SZE using either single 
dish telescopes (e.g. Mason, Myers & Readhead 2001; 
Pointecouteau et al. 2001 & 2002; De Petris et al. 2002 or 
interferometric techniques (e.g. Jones et al. 2001; Reese et 
al. 2001 & 2002; Udomprasert, Mason & Readhead 2001) 
have improved dramatically in recent years (see Birkin- 
shaw 1999 for a recent review). Measurements of the SZE 
have been used for a variety of scientific applications, in- 
cluding the estimation of the Hubble Constant (e.g. Ma- 
son, Myers, & Readhead 2001; Jones et al. 2001; Reese 
et al. 2002); determinations of the fraction, by mass, of 
baryons in the Universe (e.g. Grego et al. 2001); X-ray 
independent measurements of cluster temperatures (e.g. 
Pointecouteau et al. 2002); and constraints on cluster pe- 
culiar velocities (Holzapfel et al. 1997b, LaRoque et al. 
2002). 

The technique of interferometric SZE imaging is now 
well established. By contrast, only recently has it become 
possible to generate direct (rather than aperture synthe- 
sis) SZE images using single dish telescopes. To date, only 
two clusters, RX J1347-1145 & RX J2228+2037, both at 
z ~ 0.4, have published direct SZE images (Pointecouteau 
et al. 1999, 2001 & 2002; Komatsu et al. 2001). We 



report here a third direct SZE image. This 40 GHz (7.5 
mm) image of the (z = 0.055, Sodre et al. 1992) Abell 
cluster (Abell, Corwin, & Olowin 1989) A3667 was made 
using the Viper telescope at the South Pole as part of 
the Viper Sunyaev-Zel'dovich Survey (VSZS). The VSZS 
aims to study a complete sample of southern clusters at ra- 
dio/microwave, optical and X-ray wavelengths and it goals 
include the measurement of the Hubble Constant. We 
chose to study A3667 as our initial VSZS target because, 
as one of the brightest X-ray (io.5-2.okeV = 4.1 x 10 44 ergs 
s" 1 , David et al. 1999) clusters in the REFLEX catalog 
(Bohringer et al. 2001), A3667 is expected to have a strong 
SZE signal. A3667 is also one of the best studied clusters 
in the sky. Supporting data at other wavelengths include: 
X-ray imaging and spectroscopy (ROSAT: Knopp, Henry 
& Briel 1996; ASCA: Markevitch, Sarazin, Vikhlinin 1999 
White et al. 2000; BeppoSAX: Fusco-Femiano et al. 2001 
Chandra: Vikhlinin, Markevitch & Murray 2001 a& b 
XMM), optical data in the form of images, weak lensing 
mass maps and multi-object spectroscopy (e.g. Joffrc et 
al. 2000; Katgert et al. 1998), and radio maps (Rottgering 
et al. 1997; Hunstead et al. 1999). In addition, A3667 has 
been the focus of three-dimensional MDH/N-body numer- 
ical simulations (Roettiger, Burns & Stone 1999). 

An outline of the paper is as follows. In section 2 we re- 
view our observing strategy and data reduction methods. 
We also present the Corona map of the area around A3667. 
In section 3 we describe a joint fit, to an isothermal, tri- 
axial ellipsoidal /3-model, to the ROSAT Position Sensi- 
tive Proportional Counter (PSPC) and Corona images of 
A3667 and present an estimate of the Hubble Constant. In 
section 4, we present conclusions and discuss future plans 
for cluster observations with Viper. 
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2. THE CORONA OBSERVATION OF A3667 

The observations presented herein were made with the 
2.15 meter Viper telescope during the austral winter of 
1999. At the time, Viper was fitted with a two pixel re- 
ceiver instrument known as Corona that measured the to- 
tal incident power in the frequency range 38 to 44 GHz 
(v = 40, Coble ct al. 1999) using HEMT (high elec- 
tron mobility transistor) amplifiers. The following features 
make Viper ideally suited to the observation of low red- 
shift clusters such as A3667; its location at the South Pole, 
the two story conical shield and the flat chopping tertiary 
mirror that sweeps the telescope beam backward and for- 
ward by up to 3.6° in co-elevation. The South Pole is one 
of the driest sites in the world (Peterson et al. 2002), the 
shield minimizes ground pick up and the chopper permits 
observations of clusters that are extended over large areas 
of the sky. 

Using Corona on Viper, we have generated a 40 GHz 
image of a 2° x 3.6° region centered on A3667. We did 
so using a raster scan technique that involved taking mea- 
surements during eight complete chopper cycles at each of 
25 fixed elevation positions separated by 300" (roughly one 
third of a beam). For the A3667 observation, we used a 
3.6° chopper throw (the maximum possible) and read out 
the instrument 512 times per 2.15 s _1 chopper cycle, re- 
sulting in measurements separated by 50.6" in the azimuth 
direction. Including overheads, each raster scan map took 
a little over a minute to complete. In total, 3,103 raster 
scans were made over six contiguous days: May 28, 29, 
30, 31 & June 1, 2 1999. Every two hours during that pe- 
riod, the Carina nebula was observed as a pointing check. 
During our observation, the Corona beam pattern was well 
described by a radially symmetric Gaussian with a FWHM 
of 0.261° (15.7') and the total calibration uncertainty was 
8% (Peterson et al 2000). 

The raster scan data were initially stored as a time 
stream of detector voltages. The first stage of the anal- 
ysis was to apply a filter to the time stream. This cor- 
rected for the detector response function; removed high 
frequency noise; and eliminated any signal coupled to the 
power mains, by removing harmonics of 60 Hz. Next noisy 
and/or poor quality scans were flagged and cut from the 
analysis: Those with any > 5ct voltage spikes, those with 
tracking or chopping position errors greater than two arc- 
minutes, those for which the amplifiers were warmer than 
25K, those with excessive noise at harmonics of 60 Hz, 
those with rapidly varying voltages and those with an un- 
usually high average voltage. We note that the first data 
cut removed electrical interference in the data, and the 
latter two helped remove time periods when clouds were 
moving across the field of view. Once all seven cuts had 
been applied, 1518 scans of the initial 3103 remained or 
~36 hours of integration. 

As it is conventional to express CMB intensity changes 
as changes in brightness temperature, AT, the next stage 
of the analysis involved converting the voltage time stream 
into a AT time stream using the detector calibration. 
Then the 1518 sets of time stream were converted into 
grid maps consisting of 512 rows and 25 columns each. 
Next, the grid maps needed to be corrected for contami- 
nating signals known collectively as "chopper synchronous 
offsets" . These offsets can include signals originating from 



the ground, the atmosphere and the telescope components. 
The offset level varies with chopper position because, as 
the chopper moves, the focal plane is illuminated by differ- 
ent areas of the telescope optics (see Peterson et al. 2000). 
Under the assumption that the offsets change linearly with 
elevation, and maintained the same general waveform in 
the azimuth (chopping) direction, the offsets were removed 
from the inner regions of the scans as follows. The linear 
component in the elevation direction was removed by fit- 
ting to the top and bottom 4 rows in each of the columns. 
Next, the linear component in the azimuth direction was 
removed by fitting to the first and last 50 columns in each 
of the rows. After offset correction, the 512 x 25 points in 
the offset subtracted grid maps were assigned sky coordi- 
nates using the Peterson et al. (2000) telescope pointing 
model. Finally, the maps were converted into fully sam- 
pled images with 113" x 113" pixels. Linear interpolation 
in the elevation direction was used to create these images 
because the grid points in the raster scan were separated 
by 300" in elevation (compared to 51" in azimuth). 

In Figures Bl and B2 we present the result of averaging 
the 1,518 fully sampled images together. In these Figures, 
the data have been smoothed with a 10.3' FWHM Gaus- 
sian kernel, which gives them an effective resolution of 18' 
FWHM (the beam size of the instrument was 15.7'). The 
effective beam FWHM has been represented by the red cir- 
cle in the bottom right corner of Figure B2. In Figure Bl 
we present the full dimensions of the co-added image, in- 
cluding the edge regions used for the offset subtraction. 
The elongated blue (cold) region in the center of the image 
corresponds to the position of A3667. Other cold struc- 
tures in the map most likely correspond to fluctuations due 
to receiver noise, primary CMB anisotropies, or artifacts 
of the offset removal. It is unlikely, given the sensitivity 
and spatial resolution, that any of the cold spots in the 
map correspond to serendipitous cluster detections. Over 
the inner region of the smoothed image, i.e. the region 
not used for offset subtraction, the standard deviation per 
113" x 113" pixel was found to be ~ 60^ K. When deriv- 
ing this random noise estimate, we assumed that the noise 
was Gaussian and that each scan was independent. We 
tested the latter assumption by re-calculating the noise 
after averaging into groups of first five and then twenty 
consecutive scans. Grouping the scans in this way before 
calculating the standard deviation in each pixel did not 
change the noise estimate. This supports our assumption 
that the noise was uncorrelated in time; if it had been, 
this test would have resulted in an increase in the mea- 
sured standard deviation. 

In Figure B2 we compare the inner 1° x 1° of the Corona 
map with X-ray and radio images of the same region. Note 
that for these comparisons, the Corona data were rebinncd 
to the PSPC pixel scale (15" x 15") and that this smooths 
out the (113" x 113") pixelization apparent in Figure Bl. 
The white contours in the Figure correspond to the 0.4-2.0 
keV X-ray count rate as measured by the PSPC instru- 
ment on the ROSAT satellite. The X-ray count rate map 
was generated from data in the ROSAT public archive us- 
ing the data reduction pipeline designed for the SHARC 
Survey (Romcr ct al. 2000). This data reduction pipeline 
was based on the ESAS routines of Snowden et al. (1994) 
and includes the production of a complementary 0.4-2.0 
keV exposure map that was used during the model fit- 
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ting described below (Section 3.3). The black contours 
correspond to the interfcromctric radio map at 843 MHz 
created with the MOST interferometer (Hunstead et al. 
1999; Richard Hunstead, private communication). 

From Figure B2, it is clear that there is a significant ex- 
tended cold area close to the center of the X-ray emission. 
At the centroid (20 h 12 m 28.9 s , -56°49'51" J2000) of the X- 
ray emission, the temperature decrement in the smoothed 
Corona map is ATcmb = — 154yuK. There is also a some 
evidence for correlations between hot regions in the im- 
age and the MOST radio contours. Taken together, the 
various correlations between the Corona map and the X- 
ray and radio maps support the claim in Peterson ct al. 
(2000) that the Corona pointing reconstruction is accu- 
rate to within ±4'. The extended radio emission to the 
north-west of the cluster has been discussed in Rottgering 
et al. (1997) and Roettiger, Burns and Stone (1999). Both 
groups suggest that a merging event has produced turbu- 
lence in the outer region of the cluster which in turn ampli- 
fies magncitc fields, rcaccelerates relativistic particles and 
thus produces radio emission. Radio contamination can 
be a serious concern when creating SZE images, e.g. RX 
J1347-1145 (Pointecouteau et al. 1999, 2001; Komatsu et 
al. 2001). However, as there is no evidence of radio con- 
tamination in the region of the cluster used for the Hq 
measurement (Section 3), i.e. the core of the cluster, we 
do not discuss the radio emission from A3667 again in this 
paper. 

As mentioned above, the random noise per pixel in Fig- 
ure Bl is ~ 60/i K. In addition, the images are contami- 
nated by primary CMB anisotropies. We have estimated 
the level of the CMB contamination in the Corona images 
to be ~ 38/iK RMS using the Netterfield et al. (2001) 
anisotropy power spectrum and assuming a top hat win- 
dow function with 10 < I < 400. Therefore, some of the 
features visible in this Figure could correspond to > 2a 
CMB fluctuations. For this reason, we cannot rule out 
that the signal at the position of A3667 is free from CMB 
contamination. The northern extension of the central cold 
region could be due to such contamination. Recent obser- 
vations of A3667 by the multi-frequency Arcminute Cos- 
mology Bolometer Array Receiver (ACBAR, Runyan et 
al. 2002) on Viper should allow us to quantify the level of 
CMB contamination in the Corona map because ACBAR 
has detectors sensitive to frequencies above and below the 
null 2 in the thermal SZE spectrum. 

3. DERIVATION OF THE HUBBLE CONSTANT 

It has been known for some time (e.g. Silk & White 
1978; Cavaliere, Danese & Zotti 1979) that a measurement 
of the thermal SZE toward a cluster could result in a mea- 
surement of the Hubble Constant (Hq) when it is combined 
with X-ray observations. This is because the the SZE is 
proportional to the line of sight integral of T e n e , whereas 
the X-ray surface brightness is proportional to the line of 

sight integral of Ten 2 . These differing dependencies on 
the electron temperature and number density, allow one 
to solve for Hq, assuming we can find an accurate model 
for the distribution of these parameters along the line of 
sight. 

2 The thermal SZE results in a temperature decrement (increment) 
is no net change in CMB brightness temperature. 



We have used the A3667 data described above to con- 
strain Hq. We have performed a joint maximum likeli- 
hood (Section 3.3) analysis of the Corona and ROSAT 
PSPC images. During this analysis we have generating 
many thousands of pairs of simulated ROSAT and Corona 
A3667 images by making various assumptions concerning 
both the three dimensional distribution of electrons in the 
cluster and the performance and noise characteristics of 
the ROSAT PSPC and Viper Corona instruments (§ 3.1 
and 3.2). The technique follows that used by other authors 
e.g. Reese et al. (2002) and Grainge et al. (2002), but 
is the first application of the method to real space, rather 
than Fourier space, analysis. 

The assumption underlying our analysis was that the 
number density of electrons in A3667 follows an isother- 
mal (3, or King (King 1962), model. The /3-model can be 
parameterized by: 

/ / r \ 2 y 313/2 

n e = n e0 ( 1 + (^) J , (1) 

where n e o is the central electron number density; r is 
the radius from the center of the cluster; and r c is the 
core radius (Cavaliere & Fusco-Femiano 1978). The X-ray 
emission from A3667 is clearly not spherically symmetric 
(Figure B2), so we have further assumed that the cluster 
has a tri-axial ellipsoidal geometry (Zaroubi et al. 2001 
have shown that the ellipsoidal (3 model provides a good 
approximation to the three dimensional structure of sim- 
ulated clusters). We define the radius in Equation 1 as 
follows: 

r 2 = x 2 + {eyf + {e^l) 2 , (2) 

where the x, y, and I directions are measured along the 
principle axes of the ellipsoid. The axes arc defined so 
that x and y are in the plane of the sky and I lies along 
the line of sight. There are two isomorphic sets of models, 
those with e < 1 and those with e > 1. For the model 
fitting described below (Section 3.3), we used e < 1, i.e. 
the cllipticity was given by the ratio of the minor to ma- 
jor axes. With this convention, the x axis lies along the 
minor axis and the y axis along the major axis. We have 
parameterized the stretch along the / axis by e 1 '. If ( = 0, 
then = 1, and the stretch along the I axis is equal to 
that along the x and the ellipsoid is described as prolate. 
If C = 1, then e 1 * = e, and the stretch along the I axis is 
equal to that along the y axis and the ellipsoid is described 
as oblate. 

We have chosen to use the ellipsoidal (3 model with 
C = 0.5 for our analysis because, as shown by Grainger 
(2001), this particular model does not result in a signifi- 
cant bias in the measured H value when applied to an en- 
semble of clusters taken from N-body simulations. The (3 
model is physically motivated; it provides a reasonable, al- 
though not exact (see Sarazin 1988), description of the ra- 
dial dependence of the electron density in isolated spheres 
of isothermal gas in hydrostatic equilibrium. However, it 
does have several inadequacies. For example it predicts 
there will be a finite electron number density even at very 
large radii from the cluster center. It also includes an 

ow (above) v = 220 GHz. At the cross over, or null, frequency there 
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assumption of isothermality when certain clusters, includ- 
ing A3667, are clearly not isothermal (Vikhlinin, Marke- 
vitch & Murray 2001a). Moreover, many clusters, includ- 
ing A3667, deviate from spatial symmetry. For these rea- 
sons, future analyzes of SZE images would benefit from 
the use of non-parametric MDM-N-body simulations that 
have been constrained by observations rather than generic 
models. 

3.1. Simulated ROSAT PSPC Images 

For a cluster, the X-ray surface brightness is propor- 
tional to the line of sight integral of the square of the elec- 
tron number density. Assuming n e follows Equation 1, 
using infinite integration limits, and expressing surface 
brightness as a count rate, Cx, one obtains; 

/ n2\ "3/3+1/2 
C X (0 r ) = Cx O (l + £j , (3) 

where Cxo is the central surface brightness and 9 C and 
9 r are related to r c in Equation 1 and to r in Equation 2, 
respectively, as follows: 6 C = t c /Da and 8 r = ti-q/Da, 
where Da is the angular diameter distance and r; = o is 
given by Equation 2 when 1 = 0. 

Using Equations 2 and 3, we were able to construct sim- 
ulated ROSAT PSPC count rate maps of A3667. The 
simulated PSPC maps were generated in the same energy 
band (0.4-2.0 keV), with the same pixel size (15" x 15") and 
with the same aim point (20 h 12 m 31.20 s , -56°49'12.0") as 
the real observation (Section 2). However, rather than 
covering the full 512 x 512 pixel area of the ESAS map, 
the simulated map covered only 256 x 256 pixels (roughly 
1° x 1°). The eight inputs to the ROSAT PSPC simulation 
were as follows; 

• 0, 9 C and e to define the shape of ellipsoidal /3-model; 

• to define the projected inclination angle of cluster. 
This was defined so that <p — Q° when the major axis 
points toward the West and <f> = 90° when the major 
axis points toward the North; 

• ao and So to define the sky coordinate of the center 
of the simulated cluster. These parameters were free 
to vary and so the center of the simulated map did 
not necessary coincide with the geometric center of 
the map; 

• Cxo to define the 0.4-2.0 kcV ROSAT PSPC count 
rate in a 15" x 15" pixel located at (ao,<5o)- We note 
that this is the central count rate before the correc- 
tion for the instrument point spread function (PSF) 
had been applied (see below); 

• Cb to define a constant X-ray background in the 
same units as Cxo- 

Before embarking on a maximum likelihood analysis, 
the simulated images had to be corrected for the effects 
of the PSPC PSF. This correction was complicated by the 
fact that the PSF is a function of the radial distance from 
the center of the field of view. To account for the chang- 
ing PSF, we convolved each simulated image with 10 dif- 
ferent smoothing kernels. Each kernel corresponded to a 



model PSF appropriate to a different, evenly spaced, ra- 
dius from the aim point. The ten resulting convolved im- 
ages were then co-added using the weighting scheme out- 
lined in Equation B6. This method resulted in a dramatic 
computational speed-up, compared to a direct calculation 
of a convolution with a smoothly varying PSF model, be- 
cause it allowed us to use Fast Fourier Transforms (see 
Appendix B for more information). Even so, it would 
have been computationally intensive to convolve the en- 
tire 256 x 256 area of the simulated maps, so instead, we 
convolved only those pixels with a simulated count rate 
of least 2% of Cxo- The number of pixels meeting that 
criterion varied depending on the values of the input pa- 
rameters, but typically only comprised those within a ~ 60 
pixel radius of the center. This simplification was justified 
by the fact that cluster profiles vary slowly at large radii, 
meaning the shape of the PSF has little impact on the 
per pixel count rate outside the cluster core. We numeri- 
cally tested how much impact PSF convolution has on the 
outer regions of the cluster model and found it to be in- 
significant. To create a final simulated image, the center 
of the original 256 x 256 pixel simulation was replaced by 
the circular convolved region. We note that the edges of 
the convolved region were cropped before insertion into the 
larger map in order to correct for the effects of the artificial 
periodic boundary condition imposed by the FFT. 

3.2. Simulated Corona Images 

For every simulated ROSAT PSPC X-ray count rate 
image, we generate a complementary simulated Corona 
brightness temperature image. In the non-relativistic and 
low frequency limit, the thermal SZE brightness temper- 
ature can be calculated for a given line of sight using the 
following equation; 

AT 

= Y(xcoth(x/2)-4), (4) 

JCMB 

where Tomb is the average temperature of the CMB 
(2.73K, Fixcn et al. 1996); a; is a dimensionless observing 
frequency, x = fcT ^ IB and the Compton-Y" parameter, is 
given by the line of sight integral 

Y = cr T [ -^A^n e dl, (5) 

where o~t is the Thompson cross section and n e , T e and 
m e are the electron number density, temperature and mass 
respectively. 

Each simulated Corona image was generated with the 
same pixel size (113" x 113") and assuming the same ob- 
serving frequency (y = AOGHz) as the raster scan image 
described above (Section 2). The same input values for 0, 
8 C , e, a 0} S and <f) were used to generate both the Corona 
simulation and the companion ROSAT simulation. In ad- 
dition, input values for the electron temperature, T e , and 
the central electron number density, n e o. were required. 
The electron temperature enters the brightness tempera- 
ture calculations via the the Compton-Y" parameter (Equa- 
tion 5) and was assumed to be constant throughout the 
cluster. The central electron number density enters via 
the normalization of the /3-model (Equation 1) and was 
calculated using the expression given in Equation 7, see 
section 3.2.1. 
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The brightness temperature, AT, is proportional to the 
line of sight integral of the electron number density (Equa- 
tions 4 & 5). If n e follows a /3-model (Equation 1), one 
can obtain a simple expression for AT using infinite inte- 
gration limits that is similar to the one (Equation 3) we 
used to produce the simulated ROSAT maps, i.e.; 



AT = AT 1 + 



el 



-3/3+1/2 



(6) 



where AT is the central brightness temperature and 6 r , 
9 C and Da arc as defined for Equation 3. Unfortunately, 
as shown by Puy et al. (2000), using infinite integration 
limits leads to a significant error in the derived value of 
Hq. Therefore, rather than use Equation 6, we integrated 
Equation 1 numerically over the range — 14r c < r < +14r c . 
By using these limits we avoided included spurious SZE 
signal at large radii. We note that as a result of our def- 
inition of r (Equation 2), these integration limits defined 
an iso-density ellipsoid around the cluster. 

After the completion of the numerical integrations, the 
final steps needed to produce a simulated Corona were to 
convolve the simulated map with the Viper beam and then 
apply the same offset subtraction scheme used on the real 
data (Section 2). We performed the latter step because 
offset subtraction removes certain linear modes from the 
measured SZE signal and these modes must also be re- 
moved from the simulated images before the maximum 
likelihood analysis takes place. 

3.2.1. Calculating the Central Electron Number Density 

In Appendix A, we outline the derivation of the follow- 
ing expression for the central electron number density; 



\J CxoF \ 



' 47r(l + z) 4 eC/crQi4 £(3/3) 

D A 9 c n ^(3/3 - 1/2) ■ 



(7) 



This expression had to be evaluated before each simu- 
lated Corona map was created. The majority of the pa- 
rameters in the expression have been introduced earlier, 
but others; k, ft, F, ( and Da, require additional expla- 
nation: 

• k is the ratio of electrons to protons in the intra- 
cluster medium. For our analysis, we used k — 1.16, 
which is appropriate for a fully ionized plasma with 
30% solar abundances (Feldman 1992). 

• fi is the area over which the central X-ray count rate, 
Cxo, is measured (or simulated). For our analysis we 
used a square ROSAT PSPC pixel (15" x 15"); 

• F is a scaling factor that converts an X-ray count 
rate into an energy flux. The former can be mea- 
sured with the ROSAT PSPC and the latter can be 
calculated using theoretical emission models. The 
value of F varies with electron temperature (T e ), red- 
shift (z), metal abundance (Z) and hydrogen column 
density (njj). Assuming isothermality, each cluster 
has a characteristic F value, which can be calculated 
using an X-ray spectral package such as XSPEC (Ar- 
naud 1996). For our analysis of A3667 (Section 3.3), 
we allowed T e and ur to vary, but kept the red- 
shift and metal abundance fixed at z = 0.055 and 



Z = 0.3Zq respectively: the redshift of the A3667 is 
well known (more than 100 constituent galaxy red- 
shifts have been measured, Katgert et al. 1998) and 
the ROSAT count rate varies little with metal abun- 
dance (see Hughes & Birkinshaw 1998 for a discus- 
sion of how little effect Z has on measured values of 
H ). We created a look up table of F factors covering 
a range of T e and n# values using the raymond (Ray- 
mond & Smith 1976), wabs (Morrison & McCammon 
1983) and fake it commands in XSPEC, together 
with the appropriate ROSAT PSPC response func- 
tion. 

• C is the power to which we raise s to obtain the 
stretch of the beta model in the direction of the line 
of sight axis (Equation 2). It is impossible to con- 
strain both £ and H from the projected X-ray and 
SZE data alone. Following, Grainger (2001) we set 
£ to one half, which implies that the stretch along 
the line of sight is the geometric mean of the stretch 
along the other two axes. 

• Da is inversely proportion to the Hubble constant 
and is also a function of acceleration parameter, q n . 
During our analysis, H was free to vary, but go was 
fixed at qo = 0.5. We note that at the redshift of 
A3667, the choice of qo has negligible impact on our 
Ho determination. This is in contrast to higher red- 
shift SZE observations; Hughes & Birkinshaw (1998) 
claim that by z = 0.5, current uncertainties in the 
underlying cosmological model result in a 10-20% 
systematic error in Hq. 

In summary, under the assumption that the electron 
number density follows a (3 profile, the central density, n e0 , 
can be expressed as a function of the following parameters; 
Cxo, 6 C , (3, e, C, k, z, Da (hence H and qo) and F (hence 
riff, z, Z and T e ). Of these, z, Z, n, q and ( were fixed 
during the model fitting. 

3.3. Calculating the Joint Likelihood of the Simulated 

Maps 

After converting the simulated and observed count rate 
maps into units of counts per pixel using the ESAS expo- 
sure map (Section 2), we calculated the log likelihood of 
the simulated ROSAT PSPC images using the following 
expression, under the assumption of Poisson statistics 



ln£ PS Pc = -ln(A obs (i)!)+A obs (i) ln(7V sim (z))-iV sim (i), 

i 

(8) 

where AT sim (i) is number of counts in the i th pixel of the 
simulated map and N ^ s (i) is the corresponding number 
of observed counts. We carried out the summation over 
almost all pixels within a 1 x 1° square centered on the 
PSPC aim point. The pixels we excluded from the sum 
had been flagged as belonging to foreground or background 
X-ray point sources (see Appendix B). 

We also calculated the log likelihood function for the 
simulated Corona images assuming Gaussian noise using 
the following expression 
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ln£ C orona = -|(ln(27r|S|) + A'£ _1 A) (9) 

where the vector, A, is defined as A = d b s — d sim ; c^obs is 
a vector containing the measured brightness temperature 
values (AT b s ) in each of the pixels under test and d s j m is 
the corresponding vector of AT s i m . In Equation 9, £ is the 
covariance matrix that describes the correlation between 
the pixels under test. If the pixels were uncorrelated, then 
the covariance matrix would be populated by zeros except 
along the diagonal, where it would be populated by the 
variance in the corresponding pixel. In the case of the 
Corona map, the covariance matrix is more complex be- 
cause the pixels are correlated via atmospheric noise, the 
finite response time of the detector and the linear interpo- 
lation used during the map making process. To evaluate £ 
we chose 28 pixels lying within one beam (15.7' FWHM) 
of the center of the co-added, but unsmoothed, Corona 
map of A3667. We note that we decided to analyze data 
from this region, rather than the from the whole Corona 
map, because only near the center of the cluster did the 
SZE signal exceed the background (statistical plus CMB) 
noise. Also we wished to avoid regions overlapping with 
known radio sources. From the 1,518 individual raster 
images described above (Section 2), we compiled a 28 by 
1,518 matrix (D) with a column for each of the chosen 
pixels and a row containing the corresponding brightness 
temperature measurement. Assuming that for any given 
pixel, the 1,518 measurements were normally distributed, 
the covariance matrix for each row of D was estimated 
as follows; £ = D T CD, where D T is the transpose of 
D and C is the centering matrix. The latter is given by 
C = I — i(l„ljj). Here, I is the identity matrix; 1„ is a 
column vector of ones of length n (in this case n — 1518); 
and is the transpose of 1„. Since the 1,518 scans are 
independent and the errors are normally distributed, we 
can determine the covariance matrix for vector d b s , £, by 
dividing the covariance matrix for each row of D by the 
number of scans; i.e £ = £/1518. We note that the 28 
pixels chosen for the analysis represented only a subset of 
the total number of pixels within the central beam. The 
pixels were selected in a checker board fashion, to avoid 
including any highly correlated nearest neighbor pairs in 
the analysis. The checkerboard selection method ensured 
that the covariance matrix was non-singular (i.e. invert- 
ible) without throwing away too much data. 

The joint log likelihood of the two simulated maps, 
^n£j oint , is given by the sum of the individual log likeli- 
hoods, i.e. 

ln£ ioint = InCpspc + ln£ C orona + ItlCt, + lnC nH , (10) 

where ^n£pspc and InLcomna were defined above and 
ItiCt c and lnC nH are the log likelihoods of the electron 
temperature, T e , and hydrogen column density, nn, val- 
ues used to generate the simulated Corona map. The latter 
parameters were treated differently to the other free pa- 
rameters in our model because they have been measured by 
independent experiments. Under the assumption of Gaus- 
sian statistics, their log likelihood functions are given by 

lnC Te = -ln(a Te V^) - {T? - T° bs ) 2 / (2a 2 T J , (11) 



lnC nH = -ln(a nH V2^) - (n£ - n^ s ) 2 /(2a 2 nH ), (12) 

where and are the simulation input values and 
T° bs , n^ s 7 <iTe and a na are the measured values and their 
corresponding 1 sigma errors. In our analysis we used 
T° bs = 7.31 keV; a Tc = 0.17 keV (Markevitch, Vikhlinin & 
Murray 2001) and n°£ s = 4.71 x 10 20 cm" 2 : a nH = 1 x 10 20 
cm~ 2 (Dickey & Lockman 1990). 

3.3.1. Best fit model and errors on best fit parameters 

As described above, we were able to calculate a joint 
log likelihood for each pair of simulated A3667 maps using 
Equation 10. The most likely ellipsoidal isothermal (3- 
model for A3667 was that which maximized the joint like- 
lihood. In order to find the most likely model we created 
several thousand pairs of simulated images, adjusting the 
11 free input parameters (Cxo, Cb, C , (3, e, <ft, ao, So, Hq, 
T e and nn) slightly each time. The parameter adjustments 
were governed by the downhill simplex (Nelder and Mead 
1965) routine. We note that the first term (— ln(iV {, s (i)!)) 
in Equation 8 is independent of the model being tested, 
and therefore to speed up the calculations, was not evalu- 
ated during the /n£j int maximization. The best fit values 
for nine of the eleven free parameters in the fit are given 
in Table 1. The best fit values for the other two param- 
eters, nn and T e , were identical (to 3 significant figures) 
to the values determined by independent experiments (see 
above) . The central electron density corresponding to our 
best fit model (see Equation 7) is n e o = 4.05 x 10~ 3 cm~ 3 . 

After the most likely /3-model had been found, we used 
the likelihood ratio test (LRT) to define errors on the var- 
ious input parameters in the fit. The application of the 
LRT to SZE analysis has been described elsewhere (e.g. 
Hughes & Birkinshaw 1998; Reese et al. 2000), but we 
review the salient points here. The variable S = — 21n(£) 
has a minimum S m i n for the best fit model (i.e. when C is 
at a maximum). The Cash statistic (Cash 1979) is defined 
as C = S — S m in and is, therefore, equal to zero for the 
parameter set that defines the best fit model and greater 
than zero for any other set. The 68% confidence region for 
a particular input parameter can be estimated by finding 
the range of values that it can hold for which the Cash 
statistic that is less than unity. When making the Cash 
statistic calculations, we hold the parameter at a certain 
value and then carry out a likelihood maximization by 
varying the other parameters in the model. Several values 
of the parameter need to be tested before the C < 1 region 
can be determined. The traditional approach to the LRT 
involves measuring C, and hence C, at evenly spaced grid 
points in the n-dimensional parameter space. However, 
this approach is computationally prohibitive for a model 
with 11 parameters. We, therefore, devised a new LRT 
method which was both more efficient and more objective 
then the traditional method (it is more objective because 
it requires no a priori knowledge of the parameter error 
distributions). 

Our method proceeded as follows. First we did a 
rapid search to find values of the parameter under test 
that yielded Cash statistic values greater than unity after 
marginalizing over the other 10 free parameters. We ini- 
tiated this search by fixing the parameter at a point 50% 
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above its best fit value. We then performed a likelihood 
maximization over the remaining 10 free parameters. We 
repeated this for a fixed point 50% below the best fit value. 
If these ±50% points did not yield C > 1, then we repeat- 
edly doubled the distance between the tested values until 
they did. In practice, we only had to use this doubling 
technique for the Hubble Constant (i.e. all the other pa- 
rameters had la errors that were less than 50%). After 
finding two locations where C > 1, and knowing the loca- 
tion of the C = point, we estimated the locations of the 
C = 1 points using a third order beta spline fit. The Cash 
statistic was re-evaluated at the estimated C — 1 crossing 
points and the spline fit performed again. The process 
was repeated two more times. In this way we built up at 
least nine C measurements; i.e. one at, four above and 
four below the best fit value. With at least nine inputs 
to the spline fit, it was possible to accurately determine 
of the C = 1 crossing points. We performed the LRT er- 
ror analysis for nine of the eleven 3 free parameters in our 
model fit, see Figure B4 and Table 1. In general the the 
errors were roughly quadratic, as expected for Gaussian 
errors; however, in the case of H , the error distribution 
had a broad tail on the positive side. To demonstrate 
that the Cash Statistic was a smooth function of Hq, we 
tested some additional Hq values there were not selected 
by the root finding algorithm (represented by the crosses 
on Figure B4). Our best fit value for the Hubble Con- 
stant is 64^3q km s _1 Mpc -1 (68% confidence limits), see 
below for a discussion. Best fit values for the other pa- 
rameters can be found in Table 1. For completeness, we 
list here published core radius and [3 values derived from 
previous fits to the ROSAT PSPC observation of A3667: 
9 = 2.86', (3 = 0.54 (Mohr, Mathiesen & Evrard 1999); 
e c = 3.541^', = 0.67±g"il (Neumann & Arnuad (1999; 
model B); 9 = 4.29 ± 0.96', (3 = 0.589 ± 0.051 (Mason & 
Myers 2000). Despite the fact that the previous fits were 
made to average surface brightness profiles, rather than 
to the image directly, the published values are generally 
consistent with those in Table 1. 

The error on our H measurement is primarily due to 
the statistical instrument noise in the Corona map; we 
have ~ 40% error on the measurement of the central tem- 
perature decrement and H oc AT 2 . Other sources of 
random noise are less important because we have access 
to high quality X-ray data from ROSAT and Chandra. 
Our best fit Ho value compares favorably with other re- 
cent SZE measurements; H = 60ti(t\l) km s" 1 Mpc -1 
(Reese et al. 2002; systematic errors appear in paren- 
theses), H = 66±ii(±|) km s" 1 Mpc" 1 (Mason, Myers 
& Readhead 2001), H = H = 65^ km s" 1 Mpc" 1 
(Jones et al. 2001). Unlike our analysis, these H mea- 
surements were based on observations of multiple clusters 
(18, 7 and 5 clusters respectively). We note that errors 
of ~ ±30% and exceeding ±50% are quoted in Jones et 
al. (2001) sample and Mason, Myers & Readhead (2001) 
respectively on H measurements derived from single clus- 
ters, e.g. H = 67±28 km s_1 Mpc" 1 for the A2256 H 
measurement in Mason, Myers & Readhead (2001). These 
values demonstrates that, even though the errors on our 
H measurement are large, they are not atypical, espe- 
cially when compared to other single dish experiments. 



4. CONCLUSIONS 

As the first stage of the Viper Sunyaev-Zel'dovich Sur- 
vey (VSZS), we have presented a 40 GHz raster scan map 
of the region surrounding the z = 0.055 cluster A3667. 
The cluster was observed during the Antarctic winter of 
1999 using the Corona instrument (15.7' FWHM beam) on 
the Viper Telescope at the South Pole. The Corona im- 
age of A3667 is one of the first Sunyaev-Zel'dovich effect 
(SZE) images of a low redshift cluster. Currently, only the 
Viper telescope and the Cosmic Background Imager in- 
terferometer (Udomprasert, Mason & Readhead 2001) are 
able to make SZE images of low redshift clusters. One of 
the advantages of low redshift observations of the SZE is 
that complimentary data at other wavelengths, for exam- 
ple those necessary to make Hubble Constant estimates, 
are comparatively easy to obtain. The Corona image of 
A3367 is also one of the first direct (i.e. rather than in- 
terferometer) SZE images; only two other clusters have 
published direct SZE images. 

We have used the 40GHz map of A3667 in conjunction 
with a deep ROSAT PSPC (X-ray) image to make a max- 
imum likelihood fit to an isothermal tri-axial ellipsoidal 
(3 model. We have determined 68% confidence regions 
for nine of the eleven free parameters in the model using 
a modified Likelihood Ratio Test. Our analysis method 
includes new approaches to point source detection in X- 
ray images with varying point spread functions and also 
to error estimation in multi-dimensional parameter space. 
These innovations will be applied during the analysis of 
other VSZS clusters, but are also relevant to other SZE 
experiments. Our measurement of the Hubble Constant, 
Ho = 64I3Q km s _1 Mpc -1 , is in good agreement with 
previous, SZE determined, measurements. The error on 
this measurement is large, but not atypical for single dish 
experiments, and is being driven by instrument noise in 
the Corona map. Other parameters in the fit have been 
determined to much higher precision (1-3%). 

The Corona instrument was retired at the end of 2000 
and replaced with the ACBAR instrument (Runyan et al. 
2002). ACBAR offers both improved spatial resolution 
(5' FWHM beam) compared to Corona and simultaneous 
multi-wavelength (150, 220, 280 GHz) capabilities. Using 
ACBAR, in combination with data from CTIO, XMM- 
Newton and Chandra, we have continued the VSZS and 
are working toward completion of the SZE, X-ray and 
weak lensing observations of a complete, luminosity lim- 
ited, sample of ~ 10 clusters. These observations will im- 
prove our H measurement in several ways. First, long du- 
ration ACBAR observations produce more sensitive maps 
than was possible with Corona (see Kuo et al. 2002). Sec- 
ond, we can minimize the contribution of primary CMB 
anisotropics in our SZE images by combining data at the 
three ACBAR observing frequencies (Gomez et al. 2002). 
This feature becomes important when the instrument noise 
gets down to the level of the CMB noise and it is note- 
worthy that most other SZE experiments only operate at 
a single frequency; Diabolo (Pointecouteau 1999, 2001 & 
2002), BOLOCAM (Mauskopf et al. 2000b) and SuZIE 
(Holzapfel 1997 a&b) are exceptions, but these are not op- 
timized to image low redshift clusters. Third, combining 



For the other two, njj and T e , error estimates were available from independent experiments, see section 3.3. 
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Table 1 

The results of the joint maximum likelihood fit with 68% confidence limits 

Oc £ Cko Cj e <j> a (J2000) S 

2.9 ± +0.1' 0.7T0.02 33.4 ± 0.8 2.37T0.02 0.614 ±0.009 40.0±0.8° 20 h 12 m 28.9 ± 0.1 s -56°49'51 ± 2" 
1: Count rates in the 0.4-2.0 keV band in units of xl0~ 3 ' ' 



observations of several clusters drawn from a statistically 
complete sample reduces both random errors and system- 
atic errors associated with orientation bias (e.g. Grainger 
2001). Finally, the improved spatial resolution of ACBAR 
over Corona, means that we will be able to use indepen- 
dent beams to probe different lines of the sight through 
each cluster. We only used data from the central beam 
of the A3667 Corona image for the analysis above, but by 
using multiple lines of sight we can reduce the error on 
the overall measurement. Together, these improvements 
will allow us to measure the Hubble constant with signifi- 
cantly higher precision than was possible with the Corona 
observation of A3667 presented herein. 
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APPENDIX 

DERIVATION OF THE EQUATION FOR THE CENTRAL ELECTRON NUMBER DENSITY 

We detail below how we arrived at Equation 7. We begin with a definition for the the normalization, K, of the XSPEC 
Raymond- Smith model (Arnaud et al. 1996); 



K " r " 



J n e n p dV. (Al) 



4nDl 

Here n e and n p are the electron and proton number densities respectively and is the luminosity distance to the 
cluster. If we define n to be the ratio of electrons to protons substituting Dl = (1 + z) 2 Da, we can rewrite Equation Al 
as follows 

in-i4 r 

K = - ^- / n\dV. (A2) 

4ttL>^(1 + z) 4 k J e y ' 

Expressing the volume element as dV = D A Qd0i, we can rewrite Equation A2 in terms of an integral along the line of 
sight, 

* - / < A3 > 

where is a solid angle and where the distance along the line of sight is expressed as an angle (9i = I/Da)- 
If we assume that the electron number density follows an ellipsoidal /3-model, we can write down an expression for Kg, 
the Raymond-Smith model normalization at the center of the cluster as follows: at the projected cluster centroid the 
distance simplifies to r = e^l (since x = y = 0, sec Equation 2). Thus, along the line of sight intersecting the cluster 
center, the electron number density is given by 

n e = n e0 I 1+ f^-J J , (A4) 
which is equivalent to the following if we express I and r c as angles; 

/ , F Cn x2N "3/2/3 

n e = n e0 1+ V • (A5) 



Combining Equations A3 and A5, yields 



2X -3/2/3N 2 



* =i!™ / u i + r^i i i«, ( A6) 



4tt(1 + z) 4 k 

which can be integrated using the following identity 4 and substitution; 



f 

J — ( 



(l + A 2 + B 2 )- c dA = V^ F( ^,„y 2 V + £ 2 r C+1/2 (A7) 
u = ^-^d6i = %du , (A8) 



as follows; 

4 Taken from Mason 2000, Equation C.3 
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Kn 



10- 14 D A Qn 2 e0 6 c 
4tt(1 + z) 4 ks^ 

10- 14 D A nn 2 e0 6 c T(3/3 - 1/2) 



J{l + u 2 )-^du (A9) 



K ° ~ 4^(1 + z)4>< fp • (A10) 

We can determine the scaling between K, the normalization of the Raymond-Smith model, and a ROSAT PSPC count 
rate for a partcular set of input parameters (T e , Z, z and hr) using XSPEC, see Section 3.2.1. Defining, F, to be the 
scaling factor between K and an observed (or simulated) X-ray count rate (Cx) in a pixel of size (see section 3.2.1), i.e. 

F=^, (All) 

we finally arrive at Equation 7; 



r— — 4^1 + z)%C K i i4 F (3/3) 

VCxoF y — dm — vhw-i/2)- (A12) 



REMOVING POINT SOURCES FROM ROSAT PSPC IMAGES 

We have devised a method based on Mexican hat wavelet convolutions to identify point sources in X-ray images with 
variable point spread functions (PSFs). This method is able to exploit the speed of the fast Fourier transform (FFT) by 
approximating a smoothly varying PSF with n discrete smoothing kernels. This requires 0(np ln(p)) operations, where p is 
the number of pixels in the image, compared to the 0(p 2 ) operations required for the exact convolution. As n approaches 
,Jp the method becomes exact, but, in practice, n can be much less than ^fp. We have used this approximation both to 
find point sources in the ROSAT PSPC image of A3667, see below, and to create simulated PSPC images, see Section 3.1. 
Our technique could equally well be applied to any PSPC image and, after small modifications, it could also be used to 
analyze images from other X-ray detectors. The authors have applied it successfully to ROSAT HRI and XMM images. 

The method we use to identify the point sources draws on the theoretical framework set up in Damiani et al. (1997, D97 
hereafter) and has two parts: the selection of the optimum wavelets and the identification of sources. A brief summary is 
as follows. First the optimum Mexican hat wavelets to extract point sources located at a set of n radial distances from 
the center of the PSPC FOV are selected. The X-ray image under analysis is then convolved with each of these wavelets. 
Next, a linear combination of the n convolution products is generated. The resulting combined map approximates a 
map that has been generated by convolving the X-ray image with a smoothly varying kernel (where that varying kernal 
describes the optimum wavelet at every position in the image). We note that the optimum wavelets are scaled so that 
when they are convolved with an X-ray image the result is in units of point source photon counts (see D97). For the 
analysis of the PSPC observation of A3667 described above (Section 3.3), we used n = 10 different wavelets corresponding 
to 10 evenly spaced distances from the center of the PSPC FOV. 

The Mexican hat wavelet used as our smoothing kernal is defined by: 

sjx 2 +y 2 
z(x,y,a) = 

g(x,y,a)=g(z(x,y,a)) = (2 - z 2 ^ 2 / 2 , (Bl) 

where, a is the so called scale parameter. Using an optimization procedure, we can find the kernal best suited to find 
point sources at a particular radius, i.e. the value of a that maximizes the following expression: 

c(a) = £ £ g( ^ 0,J 7°' aMM ' ) , (B2) 

i 3 

where s(i,j) is the simulated point source image positioned at the center, (io, jo), of a large array and s is normalized so 
that JV ^2 j s(i,j) — 1. We carry out this optimization for a set of n distances, r^, from the center of the PSPC FOV. 
For each of the chosen distances we maximize c(a) by varying a using the Nelder-Mead downhill simplex algorithm and 
solve for at, where k e {0, 1, ...,n — 1}. 

To generate the simulated point source images, we used the empirically determined PSPC PSF model of Hasinger ct al. 
(1993). For any given point source, this model is a function of both its photon energy spectrum and its radial distance. 
For our analysis, we made the simplifying assumption that all point sources can be described by a generic absorbed AGN 
spectrum (n# = 4 x 10 20 cm~ 2 , a = —1), since AGN's comprise the majority of the ROSAT point source population 
(Mittaz et al. 1999). 

After determining the optimal wavelet for a set of distances from the center of the PSPC FOV, we scaled the wavelets 
so that the magnitude of point source shaped objects was preserved in the convolution. The scaled Mexican hat wavelet 
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that we use is for each of our k chosen distances is 

9k{i,j) = fflfeC ( fflfe ) g ( i ' J ''° fe )- ( B3 ) 
Next, using FFT's, we convolved the X-ray image under analysis with each of the n optimum wavelets: 

Pk(i,j) = ^2^29k(l-i,m- j)v(l,m), (B4) 

I m 

where v(i,j) describes the X-ray photon count rate image. The convolved maps pk are then combined together: 

n-1 

f>{hi) = h k(iJ)pk(i,j) (B5) 

k=0 

where the weights hk are defined as follows 

r = io) 2 + (j - jo) 2 

if r < rk-i 



r-r fc _i 



. if r ,_ 1 < r & r <rk 

h k (i,j) = h k (r) = { X^-r 1 - f ~ s ^ 

I if r k <r & r < r k +i 

if r > r k +i 

and (io,jo) is the center of v(i,j). By virtue of the wavelet scaling (equation B3), the units of the combined map are 
such that if a point source was centered on position then p(i,j) would be an estimate of the photon count rate from 
that point source. This count rate can then be converted into an estimate of the number of photons collected from a point 
source, by multiplication with an appropriate exposure map (E(i,j)). Therefore, to excise point sources from an X-ray 
image, we can simply apply a photon count threshold to p(i,j) x E(i,j). 

For the A3667 analysis described in section 3.3 we excluded pixels above an 65 count threshold in p(i,j) x E(i,j) from 
the maximum likelihood analysis. We also excluded pixels within a 4 pixels radius of p(i,j) x E(i,j) > 65 count regions. 
We note that this 80 count threshold was chosen to meet the specific needs of the analysis described in section 3.3. A 
different type of analysis, a different cluster observation, or a different instrument would likely have led to a different 
choice. 
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Fig. Bl. — The 3.6° X 2° area of the Corona 40 GHz map of A3667 in units of fiK. The map is a result of averaging 1518 raster scans 
of the region, for a total exposure time of ~ 36 hours. The map has been smoothed with a 10.3' FWHM Gaussian kernel, which gives it an 
effective resolution of 18' FWHM (the beam size of the instrument is 15.7'). The pixel size is 113" X 113". The large (~ 30' in length), cold 
(blue), area in the center of the map coincides with the X-ray emission from the cluster. 
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FlG. B2. — The central 1x1° region of the Corona map m Corona have been re-binned to the ROSAT PSPC pixel scale 

(15" X 15"). The black contours represent the MOST radio (843 GHz) observation of A3667 (Hunstcad, private communication). We note the 
correspondence between the extended radio halo and hot spots in the Corona map. The white contours represent the ROSAT PSPC 0.4-2.0 
keV image of A3667. The diameter of the red circle in the bottom right corner shows the effective beam (18' FWHM) of the Corona map. 
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Fig. B3. — The 0.4-2.0 keV PSPC count rate image of A3667 (right) and the same image after (right) point sources have been masked out. 
The masked regions were not included in the maximum likelihood analysis. 
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Fig. B4. — The Cash statistic as a function of: Hubble Constant (top left); core radius (top center); f} (top left); central (0.5-2.0 kcV) 
count rate (center left); background (0.5-2.0 keV) count rate (center center); ellipticity (center right); orientation angle (bottom left); central 
right ascension coordinate (bottom center); central declination coordinate (bottom right). The C = 1 crossing points define the edges of the 
68% confidence region for each parameter. The black dots indicate the parameter values chosen by our automated error analysis software; 
notice how few values need to be tested before the C = 1 crossing points were found. For Ho we added three extra test values by hand to 
demonstrate that there were not additional minima in the C(Hq) function. 



